Classical chaos with Bose-Einstein condensates in tilted optical lattices 
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A widely accepted definition of "quantum chaos" is "the behavior of a quantum system whose 
classical limit is chaotic". The dynamics of quantum-chaotic systems is nevertheless very different 
from that of their classical counterparts. A fundamental reason for that is the linearity of Schrodinger 
equation. In this paper, we study the quantum dynamics of an ultra-cold quantum degenerate gas 
in a tilted optical lattice and show that it displays features very close to classical chaos. We show 
that its phase space is organized according to the Kolmogorov-Arnold-Moser theorem. 
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The quantum dynamics of systems presenting chaotic 
behavior in the classical limit has been a widely stud- 
ied subject in recent years, boosted in particular by ex- 
periments using laser-cooled atoms E1M 0- 0- An 
atomic-scale analog of the kicked rotor is realized by plac- 
ing laser-cooled atoms in a pulsed standing laser wave, 
a system that displays a well-known signature of quan- 
tum chaos, the dynamical localization jfj, present in time- 
periodic systems, which consists in the suppression of the 
classical chaotic diffusion by quantum interferences. An- 
other well-known signature of quantum chaos, present in 
time- independent systems, is the fact that the distribu- 
tion of energy levels takes the shape of a Wigner distri- 
bution, a behavior experimentally evidenced for example 
in the energy spectrum of Rydberg atoms in intense mag- 
netic fields. This signature clearly cannot have a classical 
counterpart. 

Quite generally, signatures of quantum chaos have no 
direct relation to the corresponding classical dynamics. 
There are fundamental reasons for this. One is that 
the notion of phase-space orbit, fundamental in classical 
dynamics, cannot be easily translated in the quantum 
world, due to the uncertainty principle. The initial con- 
ditions of a quantum system do not correspond to a single 
point in the usual phase space (qi,Pi), as AqiApi > H/2. 
One can however project the state of the system in a ba- 
sis, and define a generalized phase space formed by the 
set of amplitudes; if the basis is adequately chosen, the 
initial condition is represented by a point in such gen- 
eralized phase space. Even then, sensitivity to initial 
conditions in not present, for the more fundamental rea- 
son that the Schrodinger equation is linear. Nonlinearity 
arises in classical physics because the dynamical variables 
(qi,Pi) are also parameters of the force. In the case of the 
Schrodinger equation, the dynamical variable is the wave 
function, whereas the potential depends on (qi,Pi), and 
eventually on t, but not on the wave function. 

A Bose-Einstein condensate (BEC) is a quantum ob- 
jects described by a nonlinear quantum evolution equa- 
tion, the Gross-Pitaevskii equation (GPE). GPE has 



proved to be able to describe the BEC dynamics with 
a reasonable precision on a wide range of situations 0]. 
GPE describes the condensate as whole (neglecting the 
non-condensed fraction), and includes a nonlinear term 
representing particle-particle interactions. It is valid for 
temperatures low enough compared to the critical tem- 
perature. If we consider a cigar-shaped BEC whose trans- 
verse length L is much smaller than the longitudinal 
length, Gross-Pitaevskii equation can be reduced to one- 
dimension, and reads: 



(i) 



where H is the one-particle Hamiltonian, and (To a cou- 
pling constant describing the interaction among parti- 
cles. Uq is related to the two-particle s-wave scattering 
length a s by U = 4nh 2 a s N/(L 2 M), where N is the total 
number of atoms, M the mass of an atom, and the con- 
densate wave function is normalized to unity. Contrary 
to the one-particle Schrodinger equation, this equation 
is nonlinear in the dynamical variable ip. Note that, as 
Uq oc N, the GPE nonlinearity has a collective character: 
modeling in an analogous way many-body interactions in 
an assembly of N thermal atoms would produce a much 
smaller nonlinear term. 

The dynamics of one dimension, tilted, periodic po- 
tentials has been often considered in the recent litera- 
turejboth for free atoms H H El EH, and for BECs 
El El El- A tilted optical potential 



V = V cos( — ) + Fx 



(2) 



is generated by superposing two counterpropagating laser 
waves, one of frequency lol — k^c — 2ttc/\l and the 
other of frequency + jt). This produces a stand- 

ing wave of step d — Al/2 whose nodes are accelerated 
with an acceleration cry/2. In the (non-inertial) reference 
frame in which the nodes of the standing wave are at 
rest, an inertial force F = —Mcj/2 appears, producing 
a "tilted" potential. In such potential, quantum dynam- 
ics consists in spatial oscillations over many steps of the 



2 



potential, known as Bloch oscillations, as recently exper- 
imentally observed with both individual atoms || and 
(weakly interacting - low U ) Bose-Einstein condensates 

We shall consider, in what follows, the particular form 
of the GPE given by 



dijjfx, t) 



1 d 2 
Am ox 



+ Fx + g\i/;\ 2 V, 
(3) 

where we introduced normalized units in which space is 
measured in units of the potential step d, energy in units 
of the recoil energy Er = fi, 2 fc|/2M, time in units of 
H/Er, force in units of 2Er/\l, m = tt 2 /2, and g — 
8a s N/{L 2 k L ) = (2/X l )(U /Er) is the (ID) nonlinear 
parameter (K = 1 in such units). The eigenstates of the 
linear part of the Hamiltonian 



H 



1 d 2 

2m dx 2 



Vq cos (2-kx) + Fx, 



(4) 



are the so-called Wannier-Stark (WS) states [jij. We as- 
sume that the depth of the potential V is large enough 
to support well localized WS states, and that the dy- 
namics can be described by the lowest-energy WS state 
in each well [ijj. We note by ip n the lowest state mainly 
localized in the n th -we\\. The invariance of WS states un- 
der translations by an integer number of potential steps, 
<Pn(%) = <Po{x— n ), is related to shifts of the eigenenergies 
given by E n = Eq + nF (Eq is an energy constant that 
can be eliminated), forming the so-called Wannier-Stark 
ladder of equally spaced levels 0] , that define the Bloch 
(or Bohr) frequency of the system lob = Fd/H. Fig. 
displays an example of evolution of a BEC in a tilted po- 
tential, and the presence of an erratic behavior suggests 
the existence of chaos in a classical sense. 

In order to get a simpler description and a better un- 
derstanding of the BEC dynamics, we now introduce a 
model obtained by projecting the wave function ip over 
the WS states: 



ip(x,t) = ^c n (t)(p n (x) 



(5) 



where the c n (t) = \[F l e ldn are complex amplitudes. This 
development is justified provided the g < 1 [18], which 
corresponds, taking typical values for the parameters, to 
N < 7 x 10 4 . Reporting Eq. © in Eq. © gives the 
evolution equation: 



nFc n + g 2J Xklm c k °l °m- 
klm 



(6) 



where the Xkim are defined (choosing phases such that 
(fin is a real function) by: 



X kit) 



tfk (x)ifl (x)if m (x)if n (x)dx. 



(7) 




1000 2000 

Time (h/E R ) 



Figure 1: Erratic evolution of the populations. The BEC is 
initially prepared on twelve neighbor WS states, and the figure 
displays the evolution of the populations /„ = |(i£>„|V>)| 2 of the 
three central (most populated) wells. The result, obtained 
by numerical integration of the GPE, clearly show an erratic 
behavior. Parameters are Vo = 5, F = 0.25, and g = 0.25. 



Due to invariance under discrete spatial translations of 
the WS states one easily shows that Xk~n i-n m-n ~ 
Xkim = Xkim, and thus: 



nFc n + g S ^ j 

klm 



As WS are essentially localized in a well, we can keep only 
couplings between nearest neighbors, that is, involving 
Xooo; Xoo±i 01 • The equation of motion then simplifies 
to: 



ic n = nFc n + gxoooCn \c n \ 2 + 

g (xoo-ic*_i + Xooic* +1 ) c 2 + 

9 (2xoo-ic«-i + 2xooic„+i) |c„| 2 + 

ff(xooi |c„_i| 2 c„_i +X00-1 |c„+i| 2 c I1+ i) (8) 



where we used the invariance of the Xkim coefficients un- 
der any permutation of indexes. Fig. displays the evo- 
lution obtained from the model, for the same parameters 
and the same initial conditions as in Fig. and one ob- 
serves the same dynamical behavior. The model is thus 
able to reproduce the dynamics of the GPE. 

It is useful to write evolution equations for the popu- 



lation /„ 



and the phase 9 r , 
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Figure 2: Same as Fig. Q but the calculation was made using 
the discrete model of Eq. Jgjl. The dynamical behavior is the 
same as in Fig. Q 



In = 



2s<T g V In yln+l (In+1 + &I n ) SOI (9 n +l - &n) + 
y/ln-1 {In + 01-n-l) sin (0 n—1 @n 

) (9) 



-nF - a g I n - 
e<J, 



\Jl n +\ (In + l + 3/3 I n ) COS (6»„+i - n ) - 
\Jln-\ (3/„ + (3I n —l) COS (0 

n — 1 @n 

)\ (io) 

with (J g = g/\g\, = Xooi/xoo-i ~ -1, and we have 
rescaled the time variable as i -> t/(\g\xooo) and the 
force as F -> F/(\g\ xooo)- Parameter e = Xoo-i/xooo 
is a small since Xoo-i/xooo <C 1 B If we neglect inter- 
well particle interactions, i.e. for e = 0, all /„ = I n (0) 
are constants of motion, and n — n (O) + oj n t increases 
linearly with time with frequency uj n = —nF — o g I n . The 
variables (I n ,0 n ) have then an "action-angle" structure 
and the system is integrable. The trajectories lay on tori 
defined by the values of the constants of motion 0|. In 
usual units the frequencies uj n are 



. C^oXOOO T 
nuj B H In, 



(11) 



where the second term on the right-hand side repre- 
sents a "frequency pulling" effect due to the nonlinear- 
ity. Moreover, one can verify that Eqs. (pj and 11 Oil 
can be obtained by the usual prescription, /„ = and 

9 n = — *gY~ from the Hamiltonian given by 



H = J2 [ nFl » + y 1 " 



2a g ey/l n I n+ i (Pl n + I n+ i) cos (0 n+1 
H Q (I) +eJ2lV n (I) cos (0 n+1 - d n )] . 



-On) 



The Hamiltonian is the sum of an integrable part and an 
non-integrable "perturbation". This form strongly evokes 
the Kolmogorov-Arnold-Moser (KAM) theorem ^^|. If 
< e *C 1, the system is said to be quasi-integrable and 
the KAM theorem states that the tori are only slightly 
deformed if initially the system is far from its resonances. 
A resonance is defined by J2 n ^n^n = 0, with £ n integer. 
Note that the KAM-like structure of our system depends 
only on the smallness of the parameter e, which can be 
modified by changing the properties of the potential. 

In order to illustrate the above conclusions, consider 
the simple case in which the BEC initially projects only 
on three adjacent WS states, labeled -1,0,1. During 
the evolution, the spread of the wave function to other 
WS states is very small [2fJ. We performed a numeri- 
cal integration of the equations of motion and plotted a 
Poincarre section [of the generalized phase space (1,0)], 
corresponding to the plane (Iq, 0i — 6q), displayed in Fig. 
El which presents characteristic features of a classically 
chaotic system, despite the fact the a BEC is a quantum 
object. 
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Figure 3: Poincarre section of the generalized (quantum) 
phase space, corresponding to the plane (Jo,0i — Oo). The 
plane is defined by i_i = 0.1, 0-i = 8o, all other angles and 
actions are zero, except I\ = 0.9 — Io, imposed by the nor- 
malization condition. One observes different types of dynam- 
ics, ranging from regular to chaotic. The arrow indicates the 
1:1 resonance, and one can see around it various frequency- 
locking signatures evidenced by the presence stability islands. 
Parameters are the same as in Fig. Q 

A principal resonance is observed when the ratio of two 
frequencies is a rational number, e. g. to n /uj m — a/b (a 
and b integers). The simplest case, a = b = 1, called 
"1:1" resonance, corresponds to the resonance condition 
[cf. Eq.CEHl] 



huj B = t/oXooo (Io - h) 



(12) 



for the two neighbor sites n = 0, m = 1. It has a simple 
interpretation: the energy levels corresponding to the two 



wells become quasi-degenerate due to the nonlinear shift, 
thus favoring population exchanges. Higher order reso- 
nances affect smaller regions of the phase space. Thus, in 
order to observe chaotic behavior, the initial wave func- 
tion shall not be too smoothly distributed on the lat- 
tice: if all /„ are almost identical, only higher order res- 
onances will show up. The arrow in Fig. indicates 
the 1:1 resonance, which is surrounded by other signa- 
tures of frequency-locking, characterized by the presence 
of stability islands. The groups of three and five stability 
islands correspond resp. to the resonances 3:1 and 5:1. 
The quantum phase space has thus a KAM-like structure. 
Far from resonances regular motion is displayed (left part 
of the Poincarre section). In this case, the nonlinear- 
ity introduces population-dependent frequencies due to 
the frequency pulling effect. The emergence of these fre- 
quencies eventually blur Bloch oscillations (which would 
correspond to vertical straight lines). 

In conclusion, we have presented in this paper a quan- 
tum system that displays a behavior very similar to clas- 
sical chaos, thanks to the presence of a nonlinearity due 
to many-body interactions. As we emphasized above, 
this is not quantum chaos, at least in the most widely 
accepted definition of the term. It is closer to classi- 
cal Hamiltonian chaos, except for the fact that it is ob- 
served in a generalized phase space. The structure of the 
phase space can be interpreted in the general frame of 
the KAM theorem. This leaves open many stimulating 
questions: Are the conditions leading to Gross-Pitaevskii 
approach essential for observing "classical" chaos with 
a BEC, or the same phenomenon can manifest itself in 
less restrictive conditions? What are the necessary and 
sufficient conditions for observing "classical" chaos with 
BECs? Can one find experimentally realizable situations 
in which a quantum effect (e.g. quantum interference) co- 
exists with chaos, and how chaos will affect it? How "clas- 
sical" chaos with BECs is affected by decoherence? These 
few questions suggest that the present work opens an in- 
teresting way for investigations of the quantum-classical 
limit in mesoscopic systems. 
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